      implicit none;
--      include 'statd.for';
      _include grid;
      integer NrDAT,NzDAT,NtDAT;
      Parameter(NrDAT=IA,NzDAT=4*IA,NtDAT=JMax);
--    NrDAT=IA;   -- points on Rcyl
--    NzDAT=4*IA; -- points on Z -- more points better for flattened bodies
--    NtDAT=JMax; -- points on theta
